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Abstract 



The nonperturbative effects in the quark form factor are considered in the Wilson loop 
formalism. The properties of the Wilson loops with cusp singularities are studied taking 
into account the perturbative and nonperturbative contributions, where the latter are 
considered within the framework of the instanton liquid model. For the integration path 
corresponding to this form factor - the angle with infinite sides - the explicit expression 
for the vacuum expectation value of the Wilson operator is found to leading order. The 
calculations are performed in the weak-field limit for the instanton vacuum contribution 
and compared with the one- and two-loop order results for the perturbative part. It is 
shown that the instantons produce the powerlike corrections to the perturbative result, 
which are comparable in magnitude with the perturbative part at the scale of order of 
the inverse average instanton size. It is demonstrated that the instanton contributions 
to the quark form factor are exponentiated to high orders in the small instanton density 
parameter. 



1 Introduction 

A convenient approach which may be applied to the investigation of the infrared behavior of 
QCD is the formulation of the non-Abelian gauge theory on loop space that becomes possible by 
virtue of the close correspondence between gauge and chiral fields [1]. One of the advantages 
of this method is that it allows one to consider both the perturbative and nonperturbative 
contributions. The basic object of study in such an approach is the gauge invariant vacuum 
average of the Wilson loop operator 
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belongs to the Lie algebra of the gauge group SU (N c ), while the Wilson loop operator J>q 1 9 ! dxA i x ) 
lies (for quark lines) in its fundamental representation. 

In particular, the objects such as path-ordered contour integrals of the gauge field emerge 
naturally in the study of hard collisions of hadrons (Drell-Yan process, deep inelastic scattering, 
etc.) and the hadron form factors in the eikonal approximation [2, 3]. The use of the Wilson 
integral formalism in the heavy quarks effective theory has also been discussed (see, e.g., [4]). 
In the case of hard processes, the integration contours become infinite, since they present 
the classical trajectories of the partons participating in the collision. Considering the parton 
scattering processes at large total energies and small momentum transfers it is revealed that the 
role of the resumed soft gluons exchanges becomes significant. While in QED this problem is 
solved in closed form yielding the well-known Sudakov form factor (by virtue of exponentiation 
of the leading one loop perturbative term), in QCD the situation appears to be more complicated 
due to its non-Abelian nature and the influence of original nonperturbative effects. Nevertheless, 
it has been shown in perturbation theory that in QCD processes the exponentiation takes 
place with the exponent given by a series in the strong coupling constant. Additional soft 
contributions may appear due to various nonperturbative effects, a number of which are closely 
related to the nontrivial structure of the QCD vacuum. In the case of the instanton ensemble 
modeling the QCD vacuum, the exponentiation of the single instanton contribution can also 
be proven (see Section 4). Indeed, although the QCD vacuum plays an important role in the 
high-energy collisions, the direct investigation of these effects remains a difficult task. The aim 
of the present work is to evaluate them within the framework of the instanton liquid model of 
QCD vacuum for the simplest possible quantity - the quark form factor. 

It is generally believed that the nonperturbative structure of QCD can be well understood 
within the instanton liquid model [5]. Considering the QCD vacuum as an ensemble of in- 
stantons, one can describe a lot of the low-energy phenomena in strong interactions on the 
qualitative, as well as the quantitative level. The importance of the instanton induced effects 
of strong interactions is also supported by lattice simulations [5, 6]. The instanton picture is 
generally considered as a fruitful and perspective framework for hadronic physics. The role 
of the instantons in hard hadronic processes has been studied intensively, both theoretically 
and experimentally. The perspectives for an unambiguous experimental detection of instanton 
contributions are believed to be optimistic and very promising. The main problem arising in 
this context is the evaluation of instanton induced effects at a hard scale of a given process [7] . 
Hence, we are going to study the instanton induced contributions to the Wilson integral over 
contours of a various geometry and apply the obtained results to the analysis of nonperturbative 
effects in the high-energy hadron collisions. Recently a similar approach has been developed 
in the study of high energy elastic and quasielastic parton-parton scattering [8]. The vacuum 
averages of Wilson operators with finite-length contours were considered in the instanton model 
in [9], where the properties of the nonlocal quark and gluon condensates were discussed. 

In the present paper, we start the investigation of the instanton induced effects in the high- 
energy QCD processes. We propose an approach which allows one to evaluate the instanton 
contributions to the Wilson integrals made of several (in)finite lines containing specific cusp 
and/or cross singularities. For this purpose, we start with one of the simplest configurations, i. 
e., the angle with infinite sides that corresponds to the integration path for the Wilson operator 
describing the soft part of the quark Sudakov form factor [10]. Such configurations enter as 
constructive units in the Wilson contours for many high energy processes. 

This paper is organized as follows. In Section 2, we briefly describe the general renormal- 
ization properties of the Wilson loops with cusp singularities and find the explicit formula in 
leading order, applicable both to perturbative and nonperturbative fields. Then we analyze the 
UV and IR divergencies and calculate the perturbative part in an arbitrary covariant gauge 



2 



demonstrating explicitly the gauge invariance of the obtained (renormalized) result. Section 3 
is devoted to the study of the instanton induced nonperturbative contributions. The expression 
for them is found in the closed form in the weak-field limit in terms of the instanton profile 
function. It is shown that the instantons actually give powerlike corrections to the perturbative 
result. The magnitude of the instanton induced effects is found to be comparable with the 
perturbative part at a certain low momentum scale. In Section 4, we prove the exponentia- 
tion in the small instanton density parameter of the all-order single instanton contribution to 
the Sudakov form factor. Afterwards, we discuss the possible applications of the results ob- 
tained above to investigations of the physically interesting processes and make some conclusive 
remarks. 



2 Renormalization and perturbative contribution 

Let us consider in short the renormalization properties of the Wilson loop. It is well known, 
that evaluation of W{C) (1) encounters divergencies. The divergencies appearing for smooth 
contours without self-intersections can be removed by means of the convenient R operation 
[11, 12]: 

W R (C,g R ,^^\imW(C,g R ,^e) = limRW(C,g,e) , (3) 

where e is the dimensional regularization parameter, \i is the UV normalization scale. In this 
case, the R operation consists in the multiplicative renormalization of the coupling constant 
9 9r an d the field — > A^. 

In the more general case, the contour may contain a number of cusps, i. e., the points 
where the derivative is not smooth. In this situation, divergencies of another type arise which 
depend on the corresponding cusp angles 7«. Although they were considered as a shortcoming 
of the loop formulation of a gauge theory about two decades ago, these singularities (and the 
anomalous dimensions related to them), are shown to play an important role in the partonic 
hard processes controlling the asymptotics of the scattering amplitudes [13]. In order to remove 
these divergencies, the generalized operation has been proposed [12]. For a contour with 
one cusp C 7 this operation is reduced to 

Wr^, g R , li, Cg = lim K 1 W{C 1 , g R , li, e) , (4) 

where the contour C 7 is treated as a subtraction point. It has been proved that the K^R 
operation renormalizes multiplicatively any loop integral with finite number of cusps [12, 14]. 
The multiplicative renormalizability means that the operation acts as 

K 1 W{C 1 , g R , li, e) = Z cusp (g R , 7, e)W{C 1 , g R , li, e) , (5) 

where 

Z cusp = W~ 1 (C 1 ,g R ,fi,e) = 1 + ( — J Z„(7,/i,e) (6) 

n=l V 71 J 

is an additional renormalization constant. The coefficients Z n are the sum of simple poles plus 
some finite terms 

n ( 1 \ 

z n(l,£) = \— A *n(ri)) + (finite terms)„ . (7) 
k=i x - e ' 

A choice of the finite part of Z n to each order actually determines the K 1 subtraction scheme. 
In order to remove pole terms we shall use the modified minimum subtraction MS'-scheme 
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The first nontrivial cusp dependent term in the expansion of W{C) (1) in powers of g 2 for 
the angle with two infinite straight line rays [(Fig. la)] contains the contributions from both 
perturbative [(Fig. lb)] and nonperturbative [(Fig. lc)] fields 

W (1) (7)=wf ) ( 7 )+W'«(7). (8) 

Regardless of which part is considered - either perturbative, or nonperturbative - we will 
write the one-loop contribution in the following form: 



W {1 \l) = - 9 —^ f dxj dy u D, u (x-y) , (9) 

where C F = ^ * s convenient to present the gluon propagator D^z) in the form 

D, v {z) = b^Mz 2 ) - d^d^z 2 ) . (10) 

Here and in what follows, we use the dimensional regularization with n = 4 — 2e, e < in 
order to control the IR-divergent terms in the integrals. The remaining UV singularity (due to 
the infinitely small z 2 in the vicinity of the cusp) will be regularized by the corresponding UV 
cutoff. 

The trajectories of the incoming and outgoing quarks [(Fig. la)] may be parametrized as 
x = vis (0 < s < oo) , y = V2T (—00 < r < 0) . 
The angle between the vectors v\ and v 2 is given in the Minkowski space by 

cosh X = M 2 ) = ^ = i + ^ , -q 2 = ( P2 - Pi ) 2 , < 2 = i, (11) 

where the quark momenta are supposed to be on shell: p 2 = p 2 . = m 2 . The continuation to the 
Euclidean space is defined as [14, 15]: 

X^il- (12) 

By virtue of the interpretation of the graph Fig. 1 as an amplitude of the elastic scattering 
of an on- mass-shell one-dimensional fermion on a color singlet potential [11, 14], we have to 
consider the quantity 

IfW^.p^-p^), (13) 
where (for some technical details of the calculations, see the Appendix) 

W^)( 7 ) = - g 2 C F [(n - 2)rf 1 (0) 7 cot 7 + d 2 (0)} . (14) 

It follows from Eq. (13) that the integrals (9) in which both points x and y belong to the same 
side of the angle do not contribute to the quantity W^\^). Hence we have within the one loop 
accuracy 

W(i) = 1 - 47ra s C F (n - 2)/i( 7 )di(0) , (15) 

where 

Mt) = 7 c °t7 — 1 (16) 

is the universal cusp factor. We should emphasize here that the expression (15) holds for 
perturbative part as well as for the nonperturbative part depending on the value 6^(0). For the 
perturbative field, Eq. (15) reflects the explicit gauge invariance in the set of covariant gauges, 
since the gauge fixing parameter £ enters only in the function d 2 (z 2 ). 
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Let us consider first the perturbative part Wp\"y). By using the free propagator in the 
Euclidean space (see the Appendix) the IR-regularized value of di(0) can be written in the 
form 

d 1 (0;e,\)= { -^- Hdaa"^ , (17) 

Lb7T z Jo 

where A 2 is the IR regularization parameter. This integral diverges at the upper (UV) limit for 
e < 0, and hence we must regularize it. To this end, we may introduce the UV cutoff /i 2 , that 
corresponds to the following replacement in the denominator of the perturbative propagator: 
z 2 — > z 2 + /i~ 2 , and finally we get 

*<o;«.m/a) Hub (?)' ■ (18) 



Thus one obtains, with one loop accuracy, 



W P ( 1 ,a s ,p/\,e) = l + ^C F h( 1 )±(l-e)(^) , (19) 



and the cusp renormalization constant Z^} within the MS'-scheme: 

Z { X(l, e) = l- ^C F h(l) (J + Invr - l) , (20) 

which is in agreement with [16, 12]. Therefore, the perturbative part of the finite renormalized 
function reads 

2 

W P (i,ii) = KmZU p ( r v,e)W P ( rf ,,e) = 1 - ^C F /i( 7 )ln^ . (21) 
The one-loop cusp anomalous dimension which satisfies the RG equation reads 

+ Wp ^^ a M"> = -r CMSp 2 ) (7,«5(^))^p(7,/i,a 5 (/i)) , (22) 

where (3{gp) is the usual QCD /3-function f3(gR) = f^j^gn = const g p + 0(g B ), which can be 
found in the one-loop order [16, 12] 

r»"V>b,as) = --£-W<Pb,VL) = ^-C F h{i) . (23) 
ala.fi it 

Hence, we reproduce the Wilson operator value for the infinite contour with the Euclidean cusp 
parameter 7. 



3 Instanton contribution 

Let us estimate the nonperturbative contribution to W(C) in the instanton model. The instan- 
ton field is given by 

A,(x; p) = A;(x; p)^- = -R^rf^x - z )M^ ~ p), (24) 
z 9 
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where R ab is the color orientation matrix (a, b = 1,2,3), a a, s are the Pauli matrices, and (±) 
corresponds to the instanton or antiinstanton. The averaging of the Wilson operator over the 
nonperturbative vacuum is reduced to the integration over the coordinate of the instanton 
center zq, the color orientation and the instanton size p. The measure for the averaging over 
the instanton ensemble reads dl = dR d 4 zo dn(p), where dH refers to the averaging over color 
orientation, and dn(p) depends on the choice of the instanton size distribution. Taking into 
account (24), we write the Wilson integral (1) in the single instanton approximation in the form 

w I (C) = j r WTrexp(ta a r)\0} , (25) 

where the phase is 

<p a = R a V^ / dx^ (x - zo)Mx ' z o; P) • (26) 

We omit the path ordering operator V in (25) because the instanton field (24) is a hedgehog in 
color space, and so it locks the color orientation by space coordinates. Performing the averaging 
over the color rotations, making use of the known property of 't Hooft symbols, 

V puV pa = ^pp^va ^pa^vp i £pupa , (27) 

and performing the subtraction of self-energy part (13) we obtain the all-order single instanton 
contribution to the cusp-dependent part of Wilson loop (1): 

wi(l) = j d 4 z J dn(p) [cos 0(7, z , p) - cos 0(0, z , p)\ , (28) 
where the squared phase 2 = (p a (f> a may be written as 

<p 2 (l,z ,p)= [( v i v i) z o ~ (viZo)(vjZ ) 



X 



f 



!• 00 r -. roc r -. 

/ da V \{{-l) l+1 v t a - z ) 2 ; p] / da'tp ((-l)^ 1 ^^' - z ) 2 ; p] . (29) 
L J Jo L J 

Although the expressions (28, 29) give the complete formula for the all-order single instanton 
contribution, in what follows we restrict ourselves to the investigation of the weak field limit. 
Phenomenologically, it is assumed that the instanton size distribution is sharply peaked at a 
certain finite value. As we will see soon, this finite instanton size p provides an UV cutoff 
for the nonperturbative gluon propagator. Therefore, for the instanton field we use the same 
dimensional IR regularization with e < and the IR parameter A as in the perturbative case 3 . 
So, in this limit the leading instanton induced term reads 

r7 2\n-4 „ „ „ r d n k ~ 

wf\ 1 ) = - 9 -^ r - jdn(p) J dxp] dy v Jj^y; Aftk; p)AZ(-k; p)e-*<*-y> . (30) 
By using the Fourier transform of the instanton field 

ij(fc;p) = -y^W(fc 2 ;p), (31) 

3 Another possible approach to divergencies in the integrals over the instanton fields is to use a constrained 
instanton solution which is characterized by the profile function exponentially decreasing at large distances 
[8, 9]. 
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Eq. (30) can be written in the form of Eq. (10) with the instantonic analogue of the function 

d 1 (z 2 ): 

d^z 2 ) -> d[{z 2 ) = -4?r [dn(p)Dj(z 2 ;p,\) , (32) 



where 



Di{z 2 ; p, A) = A 



4-71 



9 2 C F 
d n k 



-ikz 



2(p\k 2 -p) 



(33) 



Above, (f(k 2 ; p) is the Fourier transform of the instanton profile function (f(z 2 ; p) and (p'(k 2 \ p) 
is it's derivative with respect to k 2 . Now using the result (15) of the previous section, we get 
the instanton contribution in the form 



wi(r, e, A) = (n - 2)/i(7) / dn(p) D/(0; e, A, p) . 



(34) 



Consider now the renormalization of the nonperturbative part for the instanton field in the 
singular gauge, 4 where the profile function is 



(p(u;p) 



z 2 (z 2 + p 2 ) 



(35) 



By using the Fourier transform of this function one gets the IR dimensionally regularized 
Dj(0; e, A, p) in the form 

D I (0;e,\,p) = 

p ir K 2+e \ 2e f°° f°° e -p 2 (ai+/3i) 

= / delicto I dSiddo n— - . (36) 

It is clear from this expression that the instanton size plays a role of the UV cutoff, since it 
provides an exponential suppression of the integrand at large (a,/3) what corresponds to the 
large momentum. Computing this integral we have 



D I (0;e,\,p) = 



P 4 tt 2 



r(i-e) 



fl rl rl 

1 + e I / / dxdydz In 
o Jo Jo 



xz + y{l — z) 
z(l-z) 



+ 0(e 2 



(37) 



Applying the renormalization procedure as described in the previous section, we find in the 
leading order the instanton contribution to the renormalization constant 



1 



Z 7 ( 7 ,e) = 1-2/1(7) (7 jdn(p)^ + A N p 



(38) 



where the finite term is 

^ 2 3 

Ajvp = 7 £ + hi7r + — - - , 
2 

and obtain the corresponding contribution to the Wilson loop 



w?\-y, A) = 1 + tt 2 /i( 7 ) J dn{p) p 4 ln(pA) . 



(39) 



(40) 



4 We consider the singular gauge, since in the regular gauges the profile function decreases insufficiently 
rapidly at infinity, and therefore the contribution of the infinitely distant part of the contour must be taken into 
account, in contrast to the singular gauge where this contribution vanishes. 
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In order to estimate the magnitude of the instanton induced effect we consider the distri- 
bution function which has been suggested in [18] (and discussed in [9] in the framework of 
constrained instanton model) in order to describe the lattice data [6], namely the following one 

" nM = 7 °* t^r exp exp t- 2m/) ' <4i) 

where the numerical constant Cn c is determined by the number of colors 

466 e ~ 1S79Nc 

^= w _ 1)!(JVe _ 2) , ~ 0.0015, (42) 

and the string tension is accepted to be a (0.44 GeV) 2 [8, 18]. Then, using the one loop 
expression for the running coupling constant 

27r o nN c -2n f 

a s(p) = - 7rr -, : , Po = \ 1 , (43) 

(J lnpA Q cD 3 

we find the instanton contribution (40) in the form [in the distribution (41), the slow varying 
logarithmic factor due to the power of the coupling as is assumed to be constant, and taken at 
the point of the mean instanton size p] 

«,F(7,A) = l + »*(7) j {^) [-j^j \n—, (44) 

where we omitted the constant term on the right hand side since it can be easily removed by 
the finite renormalization. The expression (44) shows explicitly that the instantons yield the 
powerlike corrections to the perturbative result, which is expected from a general consideration. 
Indeed, the powerlike behavior of the nonperturbative corrections to the high energy QCD 
processes may be obtained from the renormalon analysis (see, e.g., [19]). 

It is instructive to express the result (44) in terms of the mean instanton size p and the 
instanton density n calculated directly from the distribution (41), which read respectively 

_ IW2 - 3/2) 1 

9 r(A>/2-2) v 7 ^' 1 } 

n = ^(A,/2-2) (Jl_\ 2Nc (±9EeY° (W . (46) 
2 \<X S (P)J \V2^J 

Numerically, these quantities correspond to the usual mean size po and density n of the in- 
stanton liquid model, which parametrize the deltalike approximation to the distribution (41) 
(this approach, which was proposed about two decades ago [20], now gets a strong support in 
lattice data [6]): 

dn(p) = n 5(p - p )dp , (47) 
where the parameters were estimated in [5]: 

n ~ l/m~ 4 , po ~ 1/3/ra . (48) 

To compare the instanton induced and perturbative parts, which are, respectively, 



w 7 ( 7 ) — \ — K h(-f)n 2 np 4 [lng - lnp 2 p 2 
Wp^) = 1 - ^CVM7)lng , 



(49) 
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where 



^™:! )p w o.74, 



(50) 



2 [r(/? /2-3/2)] 4 

we assume that the factorization scale p (which divides the soft and hard regions of momenta 
in the factorized quark form factor) is of order of the inverse instanton size p ~ p~ x ~ 0.6 GeV. 
Then we write the total leading order contribution to the Wilson loop expectation value in the 
form 



W( 7 ,p A) 



2tt 



7 F /i( 7 )ln(p 2 A 2 ) (l + K 



C f 



(51) 



where S = ^M=ri ~ 10 is the "classical enhancement" factor with the renormalized coupling 
constant g(p) at the energy scale p ps p~ x . The ratio of the instanton correction to the 
perturbative leading term is given by 

Instanton 



= K ^J 1 ? * 0.5 , 



1-Loop Perturb. Cf 

which is estimated using the conventional value for the packing fraction 

7r 2 np 4 pa 0.1 . 



(52) 



(53) 



One can see using the main formula (44) that the strong power suppression of the instanton 
part is partially compensated by the large factor Sq Nc+1 . This means that at the energy scale 
of order of p~ l the magnitude of the instanton induced effects is comparable to the leading 
perturbative part, and must be taken into account as well. 

Let us estimate now the two-loop perturbative contribution which may be equally important 
at the chosen low scale, where the strong coupling constant is not small enough [015(0.6 GeV) ~ 
0.5]. The two-loop result has been obtained in [14] and in the limits of small and large angles 
is given by 



(7) 



-rV/,.v,4 

2 



48 111 A 2 V8 12 + 12) A 2 

&TC F N C1 [f 8 ln% - (I - i) lug 



for small 7, 
for large 7. 



(54) 



From these expressions we get that the ratio of the instanton contribution to the perturbative 
two-loop term is approximately equal to 1.4 and varies slightly (only about 10%) with changing 
of 7 over whole allowed kinematical interval. Thus, the complete consideration of the quark 
form factor at the low momentum scale must include both the two-loop perturbative part and 
the leading order instanton one, which appear to be of the same order of magnitude. 



4 Exponentiation of perturbative and instanton correc- 
tions 

Expressions (51, 54) define the first terms of the Wilson loop expansion in gauge fields. On the 
basis of the exponentiation theorem [17] for the non-Abelian path-ordered exponentials it is 
well known that perturbative corrections to the Sudakov form factor are exponentiated to high 
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orders in the QCD coupling constant. The theorem states that the contour average Wp(C) can 
be expressed as 



W P {C) = exp 



E(^) E c n (w)F n (w) 



n=l v n ' WeW(n) 



(55) 



where summation in the exponential is over all diagrams W of the set W(n) of the two- 
particle irreducible contour averages of nth order of the perturbative expansion. The coefficients 
C n {W) oc Ci?iV™ _1 are the "maximally non-Abelian" parts of the color factor corresponding to 
the contribution coming from a diagram W to the total expression (55) in the contour gauge, 
and the factor F n (W) is the contour integral presented in the expression for W. This means 
that the essential diagrams are only those, which do not contain the lower-order contributions 
as subgraphs and, as a result, the higher-order terms are non-Abelian. 

Let us now demonstrate how the single instanton contribution is exponentiated in the small 
instanton density parameter, treating the instanton vacuum as a dilute medium [21]. The 
gauge field is taken to be the sum of individual instanton fields in the singular gauge (24, 
35), with their centers at the points z/s. In this gauge, the instanton fields fall off rapidly at 
infinity, so the instantons may be considered individually in their effect on the loop. Moreover, 
the contribution of infinitely distant parts of the contour may be neglected and only those 
instantons will influence the loop integral, which occupy regions of space-time intersecting with 
the quark trajectories. Since the parametrization of the loop integral along rays of the angle 
plays the role of the proper time, a time-ordered series of instantons arises and has an effect 
on the Wilson loop. Thus, the contribution of n instantons to the loop integral Wi{^) can be 
written in the dilute approximation as 

W$ n \-y) = Tr (u 1 U 2 ...U n U n l..U^U 1 ^) , (56) 

where the ordered line integrals U^s 

U 3 {l) = T jexp (ig da vfA^a - zj^j exp (ig J dr v^A^{v 2 t - Zj^j J 

are associated with individual instantons with the positions z/s. Because of the wide separation 
of the instantons in the dilute phase and rapid falloff of fields in the singular gauge, the upper 
and lower limits of the line integrals are extended to infinity. The line integrals C/^'s take into 
account the part of the contour that goes at infinity from +oo back to — oo and in the singular 
gauge = 1. For C/ J (C/ J '), the integral is taken over the increasing (decreasing) time piece of 
the loop. 

Then, the expression is simplified when averaging over the gauge orientations of instantons. 
The averaging is reduced to substitution of W by gjU^gJ 1 , where gj is an element of colour 
group, and independent integration of each gj over the properly normalized group measure is 
performed. Under this averaging one gets 

which is just the single instanton contribution w^^) as it is given by Eqs. (25, 28). But then 
if the averaging is done in the inverse order, from n down to 1, the entire loop integral collapses 
to a product of traces 



Wr(7)-Um n< J (7). (58) 
j'=i 
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Since the individual instantons are considered to be decoupled in the dilute medium, the total 
multiple instanton contribution to the vacuum average of the Wilson operator simply exponen- 
tiates the all-order single instanton term 1^/(7) in (28), and one has 



Thus, we have proved that in the dilute regime, the full instanton contribution to the quark 
form factor is given by the exponent of the all-order single instanton result. The exponentiation 
arises due to taking into account the many-instanton configurations effect. As it is well known, 
in QED there occurs the exponentiation of the one-loop result due to Abelian character of the 
theory. In the instanton case, the analogous result takes place since instantons belong to the 
£77(2) subgroup of the £77(3) color group and the path-ordered exponents coincide with the 
ordinary ones. Exponentiation applied to the lowest order terms (51, 54) leads to vanishing 
the Wilson loop as a power of an IR cutoff A. Then, the additional instanton contributions 
suggest that the Wilson loop decreases in the instanton medium faster with A than in purely 
perturbative theory. 



We have calculated the instanton contribution to the soft part of the quark form factor, de- 
scribed in terms of the vacuum expectation value of the Wilson loop for the contour of a special 
form (28, 29). Further analysis of high energy behavior of this form factor must take into 
account the contributions of the hard and collinear parts as well. In order to do this, we would 
need to "free" the scale \x and consider it as a factorization parameter lying between the hard 
scale Q and the IR cutoff A (A < // < Q). Then we would find that the instantonic part, being 
UV finite due to the finite instanton size, does not contain an explicit dependence on the UV 
cutoff [/,, and therefore does not contribute to the cusp anomalous dimension. However, this 
issue is out of scope of the present paper, and we did not address this question here. 

The situation may be changed essentially in some other hadronic hard processes. For exam- 
ple, the Wilson path for (in)elastic high energy quark-quark scattering consists of two infinite 
straight lines lying in two different planes, separated by the finite transverse distance [3, 8]. 
The calculation of the Wilson operator for this contour reduces formally to the evaluation of 
the path integrals analogous to the ones considered above. The only difference is that the sides 
of these contours would belong to the different planes. The transverse distance between them 
provides then an UV cutoff for the integrals, both for perturbative and instanton fields, and the 
parameter of the dimensional regularization A becomes the only (IR) cutoff for them. There- 
fore, the corresponding anomalous dimension would contain the instanton induced part, and 
the instantons would play a nontrivial role in the solution of the corresponding renormalization 
group equations. This topic will be a subject of the forthcoming study. 

Moreover, in the instanton calculations of the high energy parton scattering processes [8], 
the integrals over the instanton field can be computed in all orders by virtue of the infinite 
length of lines that form the integration path specific for this process. The latter fact leads to 
the trivial factorization of the 7-angle dependence from the corresponding Wilson integral, and 
extract the instanton contribution into the numerically calculable form factors. In our case, one 
deals with the cusp, rather than cross singularities. The all-order calculation given formally by 
Eqs.(28, 29) is quite complicated for the contour with cusp, so we postpone the solution of this 
task. 

To summarize, within the instanton vacuum model we have developed an approach which 
allows one to calculate the nonperturbative contributions to the Wilson integrals over the infi- 




(59) 
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nite contour with a cusp that represent, e. g., the classical trajectories of partons participating 
in hard collisions. We have proved that in the dilute regime, the full instanton contribution to 
the quark form factor is given by the exponentiated all-order single instanton result, see (59). 
In the weak-field limit, the instanton contribution to the soft part of the color singlet quark 
form factor is found explicitly in terms of the instanton profile function in the singular gauge. It 
is shown that the instanton induced effects are of a power type (44), but nevertheless they are 
comparable in magnitude to the perturbative ones at the scale of order of the inverse average 
size of the instanton in the instanton vacuum, see (51). 
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7 Appendix 



Below we write down some useful formulas, that have been used in the calculations. The 
dimensionally regularized free propagator in the Euclidean space reads (n = 4 — 2e): 

(A 2 vr) (2 ~ w/2) rrr:/0 (l-g/2) , 7 (n/2-l) \ 

= 4^ r W 2 " !) (V z 2(n/2-l) " z 2(n/2) J " ( la ) 

The form factors d\^(u) defined in Eq. (10) are given by (u = z 2 ) 

J . . (A 2 7r)( 2 -"/ 2 )r(n/2-2) i / \ /- j / \ 

The partial derivatives in terms of the derivative with respect to the interval u in n-dimensional 
space-time 

d 2 z = 2nd u + Aud 2 u , (2a) 
d tI d u = 25 flu d u + 4z fl z l/ d 2 u . (2b) 
Scalar products in the Euclidean space 

v l v l^u = cos7 , v^vlz^Zy = z 2 cos7 + st sin 2 7 . (3) 

For calculations of the perturbative part we used the following integrals: 

rdsdt e -^ 2 +' 2 + 2s * cos ^) = —-L (4) 
jo 2a sin7 
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and 

Hdsdt st e -(^ 2 +- cos 7 ) = I - 7cot 7 

■/o 4a 2 sin / 7 

Therefore, one has for a function d{u): 



d\ k \u) = (-) k rdaa k e- au d t (a) (6) 
Jo 

and gets 

POO ry 

^ dsdt d'(u) = -7^-^(0) , (7a) 

J q dsdt ud"(u) = 7^:^(0) , (76) 

Hdsdt std"(u) = 1 ~ 7C 9 ° t7 rf(0) . (7c) 
Jo 4sm 2 7 
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Fig. 1: The notations for the quark momenta and the total cusp- dependent part of the Wilson 
loop integral for the quark form factor (a); the leading order contributions of the perturbative 
(b) and nonperturbative (single-instanton) (c) fields; (d) the all-order single instanton result; 
(e) the exponentiation of the single instanton result. 
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